Quantum point contact as a probe of a topological superconductor 
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We calculate the conductance of a ballistic point contact to a superconducting wire, produced 
by the s-wave proximity effect in a semiconductor with spin-orbit coupling in a parallel magnetic 
field. The conductance G as a function of contact width or Fermi energy shows plateaus at half- 
integer multiples of 4e^//i if the superconductor is in a topologically nontrivial phase. In contrast, 
the plateaus are at the usual integer multiples in the topologically trivial phase. Disorder destroys 
all plateaus except the first, which remains precisely quantized, consistent with previous results 
for a tunnel contact. The advantage of a ballistic contact over a tunnel contact as a probe of the 
topological phase is the strongly reduced sensitivity to finite voltage or temperature. 

PACS numbers: 73. 23. Ad, 73.23.-b, 74.25.fc, 74.45. +c 
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I. INTRODUCTION 

Massless Dirac fermions have the special property that 
they can be confined without the energy cost from zero- 
point motion. In graphene, this manifests itself as a Lan- 
dau level at zero energy, without the usual ^hujc offset 
[l[ . The zeroth Landau level contributes half as much to 
the Hall conductance as the higher levels (because it is al- 
ready half- filled in equilibrium) , leadingto the celebrated 
half- integer quantum Hall plateaus [1,[3[. In a semiclas- 
sical description, the tt phase shift at turning points, re- 
sponsible for the zero-point energy, is canceled by the 
Berry phase of tt, characteristic for the periodic orbit of 
a Dirac fermion. 

The same absence of zero-point energy appears when 
Dirac fermions are confined by superconducting barri- 
ers, produced by the proximity effect in a topological 
insulator 0, Because of particle- hole symmetry in 
a superconductor, a state at zero excitation energy is a 
Majorana bound state, with identical creation and anni- 
hilation operators. A superconductor that supports Ma- 
jorana bound states is called topological [1, Q- 

Tunneling spectroscopy is a direct method of detec- 
tion of a topological superconductor Resonant 
tunneling into a Majorana bound state produces a con- 
ductance of 2e^ /h, while without this state the tunneling 
conductance vanishes The tunneling resonance be- 
comes broader if the tunneling probability is increased, 
and one might surmise that the resonance would vanish 
if the conductance is measured via a ballistic contact. 
We show in this paper, by means of a model calculation, 
that the contrary is true: The signature of the topolog- 
ical phase is more robust when measured by a ballistic 
contact than by a tunnel contact. 

Our model calculation is in accord with general theo- 
retical considerations [l2l . [isj , but may appear counter- 
intuitive. After all, the Majorana bound state no longer 
exists as an individual energy level if it is connected by a 
ballistic contact to a normal metal, since the level broad- 
ening then exceeds the level spacing. As we have found, 
the topological phase of the superconductor still mani- 
fests itself in the conductance of a ballistic point contact. 
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FIG. 1: Solid curves: conductance of a ballistic normal-metal- 
superconductor (NS) junction, with the superconductor in a 
topologically trivial phase (blue curve, A = 8 Ebo) or non- 
trivial phase (red curve, A = AEso)- The black dashed curve 
is for an entirely normal system (A = 0). The data is ob- 
tained from the model Hamiltonian [l^.[la| of a semiconduct- 
ing wire on a superconducting substrate in a parallel magnetic 
field (Zeeman energy Ez ~ 6i?so), for the ballistic point con- 
tact geometry shown in the inset (not to scale, d — 2.5 ho, 
W — lao)- By varying the potential Vqpc at constant Fermi 
energy Ef = 120i?ao, the point contact width w is varied 
between and W. The dotted horizontal lines indicate the 
shift from integer to half-integer conductance plateaus upon 
transition from the topologically trivial to nontrivial phase. 



in a way reminiscent of the half-integer quantum Hall 
plateaus. 



II. INTEGER VERSUS HALF-INTEGER 
CONDUCTANCE PLATEAUS 

We consider the model Hamiltonian [l3, [l^ of a 
two-dimensional semiconducting wire with an s-wave 
proximity-induced superconducting gap A. (See App. 
\X\ for a detailed description.) We have calculated the 
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scattering matrix of a quantum point contact (QPC) in 
the normal region (N) at a distance d from the supercon- 
ducting region (S), by discretizing the Hamiltonian on a 
square lattice (lattice constant a = Zso/40, with Iso the 
spin-orbit scattering length) . Our key result is presented 
in Fig. [T] The number of propagating modes in the point 
contact (and hence the transmittance Tqpc) is varied by 
changing the electrostatic potential Vqpc inside the point 
contact, at constant Fermi energy Ep. Spin degeneracy 
is removed by the Zeeman energy Ez = ^gfisB in a 
magnetic field B (parallel to the wire), so that when the 
entire system is in the normal state (A 0) the conduc- 
tance increases step wise in units of e^//i (black dashed 
curve, showing the step wise increase of the transmit- 
tance from Tqpc — 0, for a fully pinched off contact, to 
Tqpc = 8, for a maximally open contact). 

The conductance G of the NS junction is obtained from 
the Andreev reflection eigenvalues i?„ at the Fermi level. 



G — —— ^ Rn{Ep). 



(1) 



The factor of two accounts for the fact that charge is 
added to the superconductor as Cooper pairs of charge 
2e. (The spin degree of freedom is included in the sum 
over n.) The superconductor can be in a topologically 
trivial {Q = 1) or nontrivial [Q — —1) phase, depending 
on the relative magnitude of Ez, A, and the spin-orbit 
coupling energy Eso = ft^/uieslso- The blue and red solid 
curves show these two cases, where the topological quan- 
tum number Q — sign Det r was obtained in an indepen- 
dent calculation from the determinant of the reflection 
matrix As we see from Fig. [l] the conductance 

shows plateaus at values Gp, p = 0,1,2, .. ., given by 
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p if g = 1, 

P+l/2 if -1. 



(2) 



The sequence of conductance plateaus in the topolog- 
ically trivial and nontrivial phases can be understood 
from basic symmetry requirements. As discovered by 
Beri [l3| . particle- hole symmetry requires that the i?„'s 
at the Fermi level are either twofold degenerate or equal 
toOorl. (See App.lBjfor a derivation.) A nondegenerate 
unit Andreev reflection eigenvalue is therefore pinned to 
exactly this value and contributes to the conductance a 
quantized amount of 2e^//i. This is the signature of the 
topological superconductor which persists even after the 
Majorana bound state has merged with the continuum 
of states in the normal metal contact. 

If we include only the degenerate i?„'s in the sum over 
n (indicated by a prime, '^'), we may write 
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G^^(l-Q + 4y ' Rr 



(3) 



A new mode that is fully Andreev reflected thus adds 
ie^/h to the conductance, with an offset of or 2e^ /h in 
the topologically trivial or nontrivial phases. The result- 
ing conductance plateaus therefore appear at integer or 
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FIG. 2: Same as Fig. [T] but now in the presence of disorder 
(for two values of the disorder strength). The first conduc- 
tance plateau in the topologically nontrivial phase remains 
precisely quantized. 



half-integer multiples of Ae^ /h, depending on the topo- 
logical quantum number, as expressed by Eq. ^ and 
observed in the model calculation. 

The quantum point contact conductance plateaus in 
the topologically nontrivial phase occur at the same half- 
integer multiples of Ae'^/h as the quantum Hall plateaus 
in graphene, but the multiplicity of 4 has an entirely 
different origin: In graphene, the factor of four accounts 
for the twofold spin and valley degeneracy of the energy 
levels, while in the NS junction there is no degeneracy 
of the energy levels. One factor of two accounts for the 
Cooper pair charge, while the other factor of two is due 
to the Beri degeneracy of the non-unit Andreev reflection 
eigenvalues. 



III. EFFECT OF DISORDER 

While in the quantum Hall effect all plateaus are insen- 
sitive to disorder, in the NS junction this applies only to 
the first plateau. As follows from Eq. (|3]), the first plateau 
at G = (1 — Q){e^/h) is determined by the topological 
quantum number Q, which is robust against perturba- 
tions of the Hamiltonian. No such topological protection 
applies to the higher plateaus, since i?„ can take on any 
value between and 1 in the presence of disorder. 

This is demonstrated in Fig. [2J where we have added 
disorder to the model calculation (both in the normal 
and in the superconducting region), by randomly chos- 
ing the electrostatic potential at each lattice point from 
the interval [— C/disoidcr, KlisordGr]- The mean free path 
Imip [^disorder depends rather sensitively on the disor- 
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FIG. 3: Differential conductance at different values of Ef — 
Vqpc (listed in units of -Bso), for two values of the distance 
d between quantum point contact and superconductor. The 
data is taken on the first conductance plateau in the topolog- 
ically nontrivial phase (A = iEso, (/disorder = 90i5so). The 
quantum point contact is in the tunneling regime for the blue 
curve (transmittance Tqpc ~ 0.1) and in the single-mode bal- 
listic regime for the black curve (Tqpc ~ 1). The width of 
the conductance peak increases both upon increasing Tqpc 
and upon decreasing d. 



der strength. We show results for f/disordor = 90 £"30 and 
110 Eso, when the mean free path (calculated in Born ap- 
proximation) is estimated at Imfp — 9 Iso and 6 Iso, respec- 
tively. (The topologically nontrivial phase itself persists 



up to ^infp = 3lso-) 



IV. EFFECT OF FINITE VOLTAGE AND 
TEMPERATURE 

These are all results in the limit of zero applied voltage 
V and zero temperature T. There is then no qualitative 
difference between the 2e^ /h conductance resonance in 
the tunnehng regime or in the ballistic regime. A sub- 
stantial difference appears at finite voltages or tempera- 
tures. 

Considering first the effect of a nonzero applied volt- 
age, we show in Fig. [3] the differential conductance 
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J^RniEp + eV). 



(4) 



The peak centered at = is the signature of the topo- 
logically nontrivial phase The height 2e^/h of this 
peak remains the same as Tqpc is raised from to 1 
by opening up the point contact, but the peak width 
increases. For a given Tqpc, moving the point contact 
closer to the superconductor also has the effect of increas- 
ing the peak width (right panel in Fig. 




FIG. 4: The solid curve is the same data as the black curve in 
the left panel of Fig.EK^F - Vqpc = WE^o, Tqpc ^ 1), but 
on a larger voltage scale to show the resonances beyond the 
conductance peak centered at V = 0. (The curve is ±V sym- 
metric.) The dashed curve shows that the conductance peak 
becomes a conductance dip when a second mode opens up in 
the quantum point contact {Ef — Vqpc = 20 Eso, Tqpc ~ 2). 



These considerations apply to the transition from the 
tunneling regime (Tqpc 1) to the ballistic regime with 
a single transmitted mode (Tqpc « 1). If we open the 
point contact further, a second mode is partially trans- 
mitted and at Tqpc ~ 1.3 the conductance peak switches 
to a conductance dip. Fig. 0] contrasts the inverted reso- 
nances at Tqcp equal to 1 (conductance peak) and equal 
to 2 (conductance dip). The voltage scale in this figure 
is larger than Fig. [31 to show also the higher-lying reso- 
nances. 

A simple estimate for the width S ~ ?i/Tdweii of the 
conductance peak in the tunneling regime equates it to 
the inverse of the dwell time Tdwcii of an electron (effective 
mass TTicff ) in the region (of size W xd) between the point 
contact and the NS interface. For the relatively large 
mean free paths in the calculation (Zmfp > W,d), the 
dwell time for point contact widths w <^ W, d is given by 
Tdwoii — TncgWd/hTqpc, so we estimate 



m^sWd ~ Wd 



^QPC^^so 



(5) 



This formula (without numerical prefactors) qualitatively 
accounts for the increase of d with decreasing d and with 
increasing Tqpc in the tunneling regime Tqpc ^ 1, but 
for a quantitative description of the ballistic regime, in- 
cluding the switch from peak to dip, a more complete 
theory is needed. 

A similarly different robustness in the tunneling and 
ballistic regime appears if we consider the effect of a 
nonzero thermal energy fc^T on the 2e^ /h conductance 
plateau. The finite-temperature conductance is calcu- 
lated from 
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FIG. 6: Contour plot of the quantum point contact potential 
(|A3|) . for the parameters I — 0.2/so, fi^a: = IBEso, ^ujjy = 
25-Bso, Vqpc = 55i5so- This is the constriction used in the 
calculations of the conductance. 



FIG. 5: Conductance in the topologically nontrivial phase 
for different values of the thermal energy ksT. The 2e^//i 
plateau is suppressed at the smallest temperatures in the tun- 
neling regime, and only for larger temperatures in the ballistic 
regime. 



We show in Fig. [S] how raising the temperature sup- 
presses the 2e^//i conductance plateau in the topologi- 
cahy nontrivial phase. The characteristic temperature 
scale for the suppression is fc^T ~ 5, so the plateau per- 
sists longest for Tqpc ~ 1, when the line width 6 of the 
resonance is the largest. 



V. CONCLUSION 

In conclusion, we have presented a model calculation 
that shows how a quantum point contact can be used to 
distinguish the topologically trivial and nontrivial phases 
of a superconducting wire. The 2e^//i conductance reso- 
nance in the tunneling regime Q persists in the ballistic 
regime, with a greatly reduced sensitivity to finite volt- 
ages and temperatures. The characteristic temperature 
scale (for a typical value E^o = 0.1 meV of the spin-orbit 
coupling energy in InAs) reaches the 100 mK range in 
the ballistic regime, which is still quite small but within 
experimental reach. 

As more and more modes are opened in the ballis- 
tic point contact, new conductance plateaus appear at 
multiples of 4e^//i which are integer in the trivial and 
half-integer in the nontrivial phase. This sequence of 
plateaus is a striking demonstration of the role which 
the degeneracy of Andreev reflection eigenvalues plays in 
the classification of topological superconductors [H, [l^ . 
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Appendix A: Model Hamiltonian 

Our model calculations are based on the Hamiltonian 
of Refs. [13) [3 J which describes an InAs nanowire on an 
Al or Nb substrate. The Bogoliubov-De Gennes Hamil- 
tonian 



n = 
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Hr — Ep 
A* 



Ep 



A 



1 
a„ 



[Hr - Ep Ac 
[ A*<jy Ep- 



m 



(Al) 



couples electron and hole excitations near the Fermi en- 
ergy Ep through an s-wave superconducting order pa- 
rameter A. (We have made a unitary transformation to 
ensure that the condition for particle-hole symmetry has 
the form used in App. iBl ) 

The excitations are confined to a wire of width W in 
the X — y plane of the semiconductor surface inversion 
layer, where their dynamics is governed by the Rashba 
Hamiltonian 



H 



U{r) 



2moff ■ - - ' ■ h ' 2^ 

(A2) 

The spin is coupled to the momentum p = —ihd/dr 
by the Rashba effect, and polarized through the Zee- 
man effect by a magnetic field B parallel to the wire 
(in the x-direction). Characteristic length and energy 
scales are Iso = /rncsaso and Eso = mcsa^^/fi^. Typ- 
ical values in InAs are Iso = 100 nm, Eso = 0.1 meV, 
Ez = ^ffcffMs^ = 1 meV at B = 1 T. 

The electrostatic potential U = C/qpc + is the sum 
of a gate potential [/qpc and an impurity potential SU. 
The impurity potential SU{x,y) varies randomly from 
site to site on a square lattice (lattice constant a), dis- 
tributed uniformly in the interval [— C/disordcr, C^disordcr]- 

The gate potential UQPc{x,y) (see Fig. El) defines a 
saddle-shaped constriction of length 2£, containing a po- 
tential barrier of height Vqpc > 0, 



^^QPC 



saddle 



max[0, Vqpc + Usadd\c{x, y)] for \x\ > i, 
^Vqpc + ^TOcffW^y^ for I a; I < 

(A3a) 

-\m^si^l{\x\ - ef + imoffw^y^. (A3b) 
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The center (0, 0) of the constriction is placed in the nor- 
mal region at a distance d from the NS interface ed, x = d. 
The characteristic width w of the constriction at the 
Fermi energy Ep > Vqpc is defined by 



I 2{Ef - Vqpc) 



(A4) 



(This is the separation of classical turning points in the 
absence of Rashba and Zeeman effects.) 

All material parameters have the same value through- 
out the wire, except the superconducting order parame- 
ter A, which is set to zero for x < d and x > L + d. The 
length L of the superconducting region if chosen long 
enough that quasiparticle transmission through it can be 
neglected (transmission probability < 10~^). 

Using the algorithm of Ref. [201 we calculate the re- 
flection matrix r of the NS junction, which is unitary in 
the absence of transmission through the superconductor. 
Andreev reflection is described by the N x N subblock 
rhe, 



reh 

rfih 



(A5) 



The Andreev reflection eigenvalues i?„ (n = 1, 2, ... A^) 
are the eigenvalues of the Hermitian matrix product 
r^^rjjg. They are evaluated at the Fermi level for the con- 
ductance ^ or at an energy eV above the Fermi level 
for the differential conductance 



Appendix B: Beri degeneracy 

We give a self-contained derivation of the degeneracy 
of the Andreev reflection eigenvalues discovered by Beri 

El mi- 

The Hamiltonian (jAip satisfies the particle-hole sym- 
metry relation 



l\ , 1 

1 oj 1 1 



-n. 



(Bl) 



For the reflection matrix r{e) at energy e (relative to the 
Fermi level) this implies 



r(-£). 



(B2) 



At the Fermi level (e = 0) the electron and hole subblocks 
in Eq. (IA5|) are therefore related by 



^hh ^eei ^e/i ^/le 



(B3) 



Let us first assume that all i?„'s are nonzero, so that 
the matrix rhe is invertible. Unitarity = 1 requires 
that rlf^r^^ + r^hr^e = 0, hence at the Fermi level 



^ = '''ee^he " 



(B4) 



is an antisymmetric matrix. (The superscript T denotes 
the transpose.) The Hermitian matrix product 



A^A = (r 



(B5) 



has eigenvalues a„ = l/i?„ — 1, n — 1,2, .. .N. 

Let be an eigenvector of A'' A with (real, non- 
negative) eigenvalue a, so A'fA^/ = a*. Then = 
(A*)* satisfies A^A^' = ~A*AA*^* = A*{A''A'i')* = 
(aA^*)* = a^*'. The eigenvalue a is therefore twofold 
degenerate, unless ^' and ^ are linearly dependent. 

If = A* for some A, then a* = A"! A^ = 
-^*(A*)* = — jAp^f, hence a = 0. So any eigenvalue 
l/Rn — 1 ^ of A^^^ is twofold degenerate, which im- 
plies that the Andreev reflection eigenvalues i?„ 7^ 0, 1 
are twofold degenerate. 

To extend the proof to the case that some i?„'s are 
zero, we regularize the inverse and consider the matrix 



rUrherl + ^)-\ (B6) 



with e a positive infinitesimal. This matrix reduces to 
the one defined in Eq. (|B4I) if r^e is invertible and is well 
defined even if it is not. Since A^^ remains antisymmetric, 
we can follow the same steps to conclude that the nonzero 
eigenvalues of AlA^ are twofold degenerate. Evaluation 
of this matrix product using the identity 



gives the expression 



AlA^ 



a 



(B7) 



(B8) 



which has eigenvalues ae,n = (1 ^ Rn)R^{Rn + e)""*- 
These are either zero or twofold degenerate, hence we 
conclude that the i?„'s are either equal to or 1 or 
twofold degenerate. 

Notice that this Beri degeneracy is distinct from the 
familiar Kramers degeneracy (although the proof goes 
along similar lines [221). Kramers degeneracy is a con- 
sequence of an anti-unitary symmetry which squares to 
— 1. The particle-hole symmetry operation 
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ph 
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X complex conjugation (B9) 



is anti-unitary, but squares to 

In the absence of time-reversal and spin-rotation sym- 
metry, only the Beri degeneracy of the Andreev reflection 
eigenvalues is operative. This is the case for the model 
Hamiltonian (jAip considered here (with time-reversal 
symmetry broken by the Zeeman effect and spin-rotation 
symmetry broken by the Rashba effect). As worked 
out in Ref. [l^, if one or both of these symmetries are 
present, then all i?„'s are twofold degenerate — includ- 
ing those equal to or 1. The Kramers degeneracy then 
comes in the place of the Beri degeneracy, it is not an 
additional degeneracy. 
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